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In the covariant light-front quark model, we investigate the form factors of Be decays into 
D, D* , Ds, Dl ,ric, J/tp, B , B* , Bs, Bs mesons. The form factors in the spacelike region are directly 
evaluated. To extrapolate the form factors to the physical region, we fit the form factors by adopting 
a suitable three-parameter form. At the maximally recoiling point, b ^ u,d, s transition form fac- 
tors are smaller than b ^ c and c d,s form factors, while the b u, d, s, c form factors at the zero 
recoiling point are close to each other. In the fitting procedure, we find that parameters in 
and yl^''^" strongly depend on decay constants of B* and B* mesons. Fortunately, semileptonic 
and nonleptonic Be decays are not sensitive to these two form factors. We also investigate branch- 
ing fractions, polarizations of the semileptonic Be decays. Be {rie, J/tp)liy and Be — > {Ba, B*)li' 
decays have much larger branching fractions than Be [D , D* , B , B'')lu . For the three kinds of 
Be -* Vlf decays, longitudinal contributions are comparable with the transverse contributions. 
These predictions will be tested on the ongoing and forthcoming hadron colliders. 

PACS numbers: 13.20.He, 12.39.Ki 

I. INTRODUCTION 

B meson decays provide a golden place to extract magnitudes and phases of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements, which can test the origins of CP violation in and beyond the standard model (SM). There 
has been remarkable progress in the study of semileptonic and nonleptonic B meson decays. Experimentally, the two 
B factories have accumulated more than 10® BB events. Some rare decays with branching fractions of the order 
10^^ have been observed. On the theoretical side, great successes have also been achieved: apart from contributions 
proportional to the form factors, the so-called nonfactorizable diagrams and some other radiative corrections are taken 
into account. All of them make B physics suitable for the precise test of the SM and the search of new phenomena 
(See Ref. [ij for a recent review). 

Compared with B mesons, Be meson is heavier: the mass of a BcBc pair has exceeded the threshold of T{AS), thus 
Be mesons can not be produced on the B factories. But Be meson has a promising prospect on the hadron colliders. 
The Large Hadron Collider (LHC) experiment, which is scheduled to run in the very near future, will produce plenty 
of Be events. With more data accumulated in the future, the study on Be mesons will be of great importance. Be 
meson can decay not only via the b ^ q (q=u,d,s,c) transition like the lighter Bu.d,s. mesons, but also through the 
c q (q=u,d,s) transitions. The CKM matrix element in the c ^ s transition \Ves \ ~ 1 is much larger than the CKM 
matrix element \Vcb\ ~ 0.04 in 5 — ^ c transition. Although the phase space in c ^ d, s decays is smaller than that in 
b ^ c transition, the former decays provide about 70% to the decay width of Be- This results in a larger decay width 
and a much smaller lifetime for the Be meson: tb^ < ^tb- The two heavy b and c quarks can annihilate to provide 
a new kind of weak decays with sizable partial decay widths. The purely leptonic annihilation decay Be — > li' can be 
used to extract the decay constant of Be and the CKM matrix element Vcb- 

Semileptonic Be decays are much simpler than nonleptonic decays: the leptonic part can be straightforwardly 
evaluated using perturbation theory leaving only hadronic form factors. In two-body nonleptonic Be decays, most 
channels are also dominated by the Be transition form factors. Thus the Be transition form factors have already 
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received considerable theoretical interests 
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In the present work, we 



will use the light-front quark model to analyze these form factors. The light front QCD approach has some unique 
features, which are particularly suitable to describe a hadronic bound state 17[. Based on this approach, a light-front 
quark model with many advantages is developed 
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22l |. This model provides a relativistic treatment of 



the hadron and also gives a fully treatment of the hadron spin by using the so-called Melosh rotation. The light front 
wave functions, which describe the hadrons in terms of their fundamental quark and gluon degrees of freedom, are 
independent of the hadron momentum and thus are explicitly Lorentz invariant. In the covariant light-front quark 
model 22], the spurious contribution, which is dependent on the orientation of the light-front, becomes irrelevant in 



the study of decay constants and form factors and that makes the light-front quark model more selfconsistent. This 
covariant model has been successfully extended to^investigate the decay constants and form factors of the s-wave and 
p-wave mesons 
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24l |. the heavy quarkonium [2 



Our paper is organized as follows. The formalism of the covariant light-front quark model is presented in the next 
section. Numerical results for the form factors and decay rates of semileptonic Be decays are given in Section IIIII 
We also compare our predictions of form factors with those evaluated in the literature. Our conclusions are given in 
Section HVl In the Appendix A, we give the relation between the form factors defined in various studies on Be decays 
and the widely used Bauer-Stech-Wirbel (BSW) form factors [26|. In the Appendix B, we collect some specific rules 
when performing the integration. 



II. COVARIANT LIGHT-FRONT QUARK MODEL 

Be ^ P,V (P, V denotes a pseudoscalar and a vector meson, respectively) form factors induced by vector and 
axial- vector currents are defined by 

{P{P")\V,\Be{P')) - f+{q')P, + f-{q')q,, (1) 
{V{P",e"*)K\Be{P')) = e^.o.pe"*''P^q^g{q^), (2) 
{V{P",e"*)\A^\Be{P')) = -z{e;'*/(9')+e"*-P[P^a+(g2) + g^a_(g2)]}, (3) 

where P — P' + P" , q = P' — P" and the convention 60123 = 1 is adopted. The vector and axial- vector currents are 
defined as ip^^ip' and i/'7m75V''- Iti b ^ q (q — u, d, s, c) transition, and ip' denotes the q quark field and the b quark 
field, respectively; while in c ^ g' {q' = u, d, s) transition, ip and ■0' denotes the q' quark field and the c quark field. 



respectively. In the literature, the Bauer-Stech-Wirbel (BSW) [2q| form factors are more frequently used: 



{PiP")\V,\Be{P')) = ( P, - ^^^^-^q,) Ff^^iq-) + ^^^^^^,P„^^^(g^), (4) 



{V{P",e"*)\VMP')) = z £/^-/3e"*''P"/^^^^(g'), (5) 

{V{P\e"*)\A,\Be{P')) = ^|(mB.+™v)C^^''(9') T 

-2my ^ qM^.'''i<i') - A^.^'^ie)]]- (6) 

These two kinds of form factors are related to each other via: 

F^^^'iq') ^ Mq% fo^- V) = /+(g^) + . 2 /-(g'). 

F^^^(g^) = -{ms^+my)giq% ^^^(g^) = 1^^, 

niB^ + my 

^f=^(g2) = {mB^+mv)a+{q% A^^"" {q') - A^^"" {q') ^ ^ a_{q% (7) 



-P2 

FIG. 1: Feynman diagram for Be ^ P,V decay amplitudes. The X in the diagram denotes the vector or axial- vector transition 
vertex while the meson-quark-antiquark vertices are given in the text. 

with Af-'^iO) = A^-'^iO), and 

Ai^-iq') = ^^^^^Af - a^:^^ A^^^iq^). (8) 

In the covariant light-front quark model, we will work in the (7+ = frame and employ the light-front decomposition 
of the momentum P' = (P'" , P'+, P^), where P'± = P'" ± P'^, so that P'^ = P'+P'- - Pf . The incoming and 
outgoing mesons have the momenta P' = p[ + P2 and P" = + p2 and the mass M' and A/", respectively. For the 
Be transition form factors, ill' = tub^- The quark and antiquark inside the incoming (outgoing) meson have the mass 
m'/"^ and 1712 and the momenta p'/"'' and p2 , respectively. These momenta can be expressed in terms of the internal 
variables (xi^p'^J as 

P'l^ = 2;i,2P'+, p[,2± = 2^1,2^1 ±p1, (9) 
with xi + X2 — 1. Using these internal variables, one can define some useful quantities for the incoming meson: 



Xi X2 u V u i. 

where e|'^ can be interpreted as the energy of the quark or the antiquark and Mq can be viewed as the kinematic 
invariant mass of the meson system. The definition of the internal quantities for the outgoing meson is similar. To 
compute the hadronic amplitudes, we require the Feynman rules for the meson-quark-antiquark vertices (iF^^): 

iV'p = H'pj,, (11) 
iVy = iHy[-ff,- —{p[- P2)^]. (12) 

For the outgoing meson, one should use iijo^'lilo) for the relevant vertices. 

In the conventional light-front quark model, the constituent quarks are required to be on mass shell and physical 
quantities can be extracted from the plus component of the current matrix elements. However, this framework suffers 
from the problem of non-covariance because of the missing zero-mode contributions. In order to solve this problem, 
Jaus has proposed the covariant light-front approach which provides a systematical way to deal with the zero-mode 
contributions ^22^1 . Physical quantities such as decay constants and form factors can be calculated in terms of Feynman 
momentum loop integrals which are manifestly covariant. For example, the lowest order contribution to a form factor 
is depicted in Fig. [T] and the P — > P transition amplitude is given by: 



f ,4 , H'p{H'l,) pp 
'2nY J ^^N[N'{N2 ^ ' 



= I dy^^^f^S--, (13) 
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2p' [M'2 + Af 
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2Ar2 - (m'l - TO2)^ - (to" - TO2)2 + (to'i 



_ 2Af'2 + TV^ - N'i + 2iV2 + 2(to; - - (to'i - m![f 



(14) 



In practice, we use the hght-front decomposition of the loop momentum and perform the integration over the 
minus component using the contour method. If the covariant vertex functions are not singular when performing the 
integration, the transition amplitude will pick up the singularities in the anti-quark propagator. The integration then 
leads to: 
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with = — X2qA^- The explicit forms of /i'^^ and uj^^ for the pseudoscalar and vector meson are given by: 

1 
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where ip' is the light-front wave function for pseudoscalar and vector mesons. After this integration, the conventional 
light-front model is recovered but manifestly the covariance is lost as it receives additional spurious contributions 
proportional to the lightlike four vector uj — (0,2,0^). The undesired spurious contributions can be eliminated by 
the inclusion of the zero mode contribution which amounts to performing the integration in a proper way. The 
specific rules under this p~ integration are derived in Ref. 
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2311 and the relevant ones in this work are collected in 



the Appendix B. 

Using Eqs. p4|) - (fT7)) and taking the advantage of the rules in Ref. 
form factors: 
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23|, we obtain expressions for the P ^ P 



167r3 



167r3 



dx2d p 



X2N[N{' I 
X2{m'i — m'l)'^ — xi{m[ 



Xi(Af^2 
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'X2N[N'1 



— 2:1X2 M'^ ^ p'l — m[m2 + {m'{ — m2){x2m[ + Xim2) 



P'l 



A.p'±-q±y 



M"^ - X2{q^ + q ■ P) 



-{X2 - xi)M''^ + 2xiM'^ - 2{m[ - m2){m[ + to'/) 



(18) 



Similarly, the P ^ V transition amplitudes are given by: 



p,PV _ -3 
" ~ (27r)4 
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where 



qPV _ (qPV qPV 



= -2ie^,^f3{p'{'P^{m'{ - mi) +pTq^{m'i + mi - 2m2) + g«P^mi} 



+2.g^^|m2(g2 - - ^" - mi^ - m'C^) - m[{M"^ - N'{ - N2 - m'i^ - ml) 
~m'l{M'^ - N[- N2- m!^ - m\) - 2m'^m'[m^ 

+^p'n.AAmi - mi) - 2{P^,q^ + + 2q^g„)mi + 2pi^F.(mi - mi') 
+Vi^g>/(3mi - mi' - 2m2) + 2P^pi^(mi + vn![) + 2g^pi^(3mi + m'/ - 2m2) 
+ - - P,){2pi^[M'2 + M"2 - g2 _ 27V2 + 2(mi - m2)(mi' + m^)] 

- 2M'2 + 7V{ - N'i + 27V2 - (mi + mi')^ + 2(mi - m^fX 
+P^[<z2-7V{-7V{'-(mi+mi')^]}. 

The above equations give the expression for P V form factors: 
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-2m2{m[ - 'm2)]{A\ 
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A. 
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The functions Z2 and a[^\ Af\ and Z2 are given in the Appendix B. Expressions for the BSW form 

factors can be directly obtained through the simple relation given in Eq. ([7]). 



III. NUMERICAL RESULTS 

The qq meson state is described by the light-front wave function which can be obtained by solving the relativistic 
Schrodinger equation. But in fact except for some limited cases, the exact solution is not obtainable. In practice, we 
usually prefer to employ a phenomenological wave function to describe the hadronic structure. In this work, we will 
use the simple Gaussian-type wave function which has been extensively examined in the literature: 



M = .A^, (22) 

dX2 X1X2MQ 

The parameter /?', which describes the momentum distribution, is expected to be of order Aqcd- It is usually fixed 
by meson's decay constant whose analytic expression in the covariant light-front model is given in [23l |. The decay 
constant of /j/^ can be determined by the leptonic decay width 

Tee = T[J/ip e^e ) = , (23) 

3m j/^ 

where Qc = +2/3 denotes the electric charge of the charm quark. Using the measured results for the electronic width 
of J/'0 0: 



Tee = (5.55 ± 0.14 ± 0.02)keV, (24) 

we obtain /j/^ = (416 ± 5) MeV. Under the factorization assumption, the decay constant of rjc has been extracted 
by CLEO collaboration from B rjcK decays [2^: 

/^^ = (335 ± 75)MeV, (25) 

where the central value is about 20% smaller than that of J/i^- In this work, we will assume the same decay constant 
for rje as that oi J/tp. We also introduce an uncertainty of 20% to this value. Decay constants for charged pseudoscalars 
are usually derived through the purely leptonic decays: 

r(P - l^) = ^^^^f.mfmAl - (26) 



The experimental results for the decay constants of charmed mesons are averaged as [29( 1: 

fjj^ = (273 ± 10)MeV, fo = (205.8 ± 8.9)MeV. (27) 

As clearly shown in the above equation, the uncertainties for these decay constants are less than 5%. It provides a 
solid foundation for the precise study on Be transition form factors. In the heavy quark limit, the decay constant fjj* 
of a vector heavy meson D* is related to that of a pseudoscalar meson through: 



fn^^fDxJ^, (28) 
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TABLE I: Input parameters and /3' (in units of GeV) in the Gaussian- type light-front wave function (|22p . Uncertainties of 
are from the decay constants as discussed in the text. 



'mu,d 


rUs 




mt 


0.25 


0.37 


1.4 


4.8 


P'd 


P'd- 




/^Dt 




0.366lS]-°o 


o.6ool«:°l^ 


0.438iHfo 


P'b 








A t:t:t:+0.048 
U.OOO_Q Q4g 


n c;9fi+0 033 


626+°"^^ 


0.599^nS 






/3s. 




0.8141°°^^ 


fi'?2+o °°5 

U.OOZ_g QQ5 


0.890t°:»- 





where mjj and ?7i£). denotes the mass of the pseudoscalar and vector meson, respectively. That implies /d* < fo since 
Tod* > niD- In the following, we will use the same values for the decay constant of the vectors and pseudoscalars. To 
compensate the differences, we will also introduce an uncertainty of 10% to the decay constants. Decay constants for 
the bottom mesons are employed by: 

fg^ = (400 ± 40)MeV, Jb = (190 ± 20)MeV, fs, = (230 ± 20)MeV. (29) 

These values are slightly smaller than results provided by Lattice QCD [s^l: 

Jb, = (253 ± 8 ± 7)MeV, /s^ = (489 ± 4 ± 3)MeV. (30) 

Decay constants of the vector B mesons are used as: Jb- ~ (210 ± 20) MeV and /s* = (260 ± 20) MeV which 
are about 10% larger than those of the pseudoscalar B mesons. Shape parameters /3's determined from these decay 
constants, together with the constituent quark masses used in the calculation, are shown in table H] The consistent 



quark masses are close to the ones used in Ref. j23l. |24|. To estimate the uncertainties caused by these quark masses, 
we will introduce the uncertainties of 0.03 GeV and 0.1 GeV to the light quark masses and the heavy quark masses, 
respectively. The masses (in units of GeV) of hadrons are used as [27| : 

TUB, = 6.276, viD = 1.8645, tod* = 2.0067, tod, = 1.9682, mz,. = 2.112, 
TO^^ = 2.9804, TOj/^ = 3.0969, ms = 5.279, ms- = 5.325, tub, = 5.3675, 
TOfi. — niB* + ruB, — ruB- (31) 

If a light meson is emitted in exclusive nonleptonic decays, only the form factor at maximally recoiling point (q^ ~ 0) 
is required but the (/^-dependent behavior in the full > region is required in semileptonic Be decays. Because of 
the condition q^ = imposed during the course of calculation, form factors can be directly studied only at spacelike 
momentum transfer q^ ~ ^ 0; which are not relevant for the semileptonic processes. It has been proposed in 



23| to parameterize form factors as explicit functions of in the space-like region and one can analytically extend 
them to the time-like region. To shed light on the momentum dependence, we will choose the parametrization for the 
b quark decays: 

F{q^) = F(0)exp(cis + C2s'), (32) 

where s = q'^ /rn^^ and F denotes anyone of the form factors Fi, Fq and V, Aq, Ai, A2. But for c ^ u,d,s transitions, 
we find that the fitted values for the two parameters ci , C2 are not stable and thus we adopt the optional three- 
parameter form: 

F{q^) = —. (33) 

' i_ ^ + 5(^)2 ^ ' 
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In the procedure to fit the form factors A^^^ and A^''^' , we find that the shape parameters {uifitjS) strongly 
depend on the decay constants /b» and Jb*- In this case, our predictions on these two form factors are unrehable, 
thus we refrain from predicting these two form factors. Fortunately, the ambiguity of A^''^ and A^''^' will not affect 
the physical quantities in various physical decay channels. As we can see from equation ([8]), the masses of B*,B* 
mesons are very close to that of thus the second term in the right hand side is negligible. The form factor Aq, 
which is relevant for the nonleptonic Be — > B*{B*)P decays, receives small contributions from A2. Contributions 
from A2 to the Be B*{Bl)lv decays and Be B*{B*)V decays are also small which will be shown in the following. 

Our predictions of the remanent form factors are collected in table |TT] and table IIIII The first kind of uncertainties 
shown in these tables are from those in decay constants of the Be meson and the final mesons; while the second kind 
of uncertainties are from those in the constituent quark masses. Several remarks are given in order. First, from these 
two tables, we can see that the Be D,D*,Ds,D* form factors at maximally recoiling point {q^ — 0) are smaller 
than the other ones. It can be understood as follows. In Be D, D* , D^, D* transitions, the initial charm quark is 
almost at rest and its momentum is of order me] in the final state, the meson moves very fast and the charm quark 
tends to have a very large momentum of order mf,. In this transition, the overlap between the wave functions is limited 
which will produce small values for the form factors. In Be rje, J/ip transitions, the spectator charm antiquark in 
rje, J play the same role with the charm quark generated from the weak vertex. The light-front wave function of 
the charmonium is expected to have a maximum aX E — — « me- The overlap between the initial and 
final states' light-front wave functions in Be — > ?/c, J li^ becomes larger, which certainly induces larger form factors. It 
is also similar for the Be ^ B, Bs form factors. Secondly, the Be Ds,ric form factors at the zero recoiling point are 
close to each other. The initial charm quark is almost at rest and its momentum is of order me- In these two kinds 
of transitions, the charm spectator in the final states tends to posses a momentum of order rrie. The overlaps of the 
wave functions in Be Da.rje transitions are expected to be in similar size. Thirdly, the SU(3) symmetry breaking 
effects in Be D, Ds and Be D* , D* form factors are quite large, as the decay constant of Ds is about one third 
larger than that of the D meson. But in Be B, Bg and Be ^ B*, B* transitions, the SU(3) breaking effect is small, 
because the decay constants /^(,) are in similar size. Fourthly, since the uncertainties from decay constants of 
Z?, Dg, J /i}} are very small, the relevant uncertainties to the form factors are also very small. 

In the literature, there already exist lots of studies on Be transition form factors 0, [EQBBSBfl 
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13l Il4l . llSl . |l6[ and their results are collected in table IIVI and table |Vl Since J/V' can be easily reconstructed 



by a lepton pair on the hadron collider, the Be —f J j"^ form factors have been widely studied in many theoretical 
frameworks. In a very recent paper 15], the authors have derived two kinds of wave functions for the charmonium 
state under harmonic oscillator potential and Coulomb potential. They also used these wave functions to investigate 
the Be r]e, J/ip form factors under the perturbative QCD approach. Compared with their results, our predictions 
are typically smaller. The main reason is that they have used a much larger decay constant /b^. Regardless of this 
effect, our results are consistent with theirs. Results collected in IIVI (including ours) have large differences which can 
be discriminated by the future LHC experiments. The Be Ds,D* is described as the FCNC 6 — > s transition at 
the quark level which is purely loop effects in the SM. As a consequence, this transition has a very small Wilson 
coefficient and the Be — > Dg, D* form factors are less studied in the literature. Similar with the b u, s, c transitions, 
predictions of the c s transition form factors have large differences between different methods. As indicated 

from these two tjibles, results evaluated in Refs. 0, 0, [l2, 14 1 are different with the other ones and ours to a large 
extent. In Ref. 9], all of the results except for the Be to charmonium transitions are larger than the other results: 
the authors have taken into account the as/v corrections and the form factors are enhanced by three times due to the 
Coulomb renormalization of quark-meson vertex for the heavy quarkonium Be- Moreover, small decay constants for 
the B meson are adopted which also give large form factors: fs = 140— 170 MeV, fs'/fs = 1-11 and fsJfB = 1-16. 
In Ref. [1^, the authors have chosen the chiral correlation functions to derive the form factors in the light-cone sum 
rules. Although only the twist-2 distribution amplitudes (DAs) contribute and contributions from the twist-3 DAs 
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TABLE II; Be D, D* , Ds,Ds, rjc, J/V* form factors in the light-front quark model. The uncertainties are from the Be decay 
constants and the decay constant of the final state mesons. 
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08+° 
— u 


01+0 
01-0 


01 
01 





421° 


02+0 
01-0 


05 
04 


3 


18lo 


28+0 
25-0 


24 
23 





65+2 
— u 


06+0 
04-0 


06 
05 




D' 





071° 


01+0 
01-0 


01 
01 





5llg 


01+0 
01-0 


07 
06 


3 


781° 


26+0 
24-0 


23 
23 





8012 


04+0 
04-0 


06 
05 


pB, 


Ds 





281° 


02+0 
02-0 


02 
02 


1 


24lg 


04+0 
05-0 


09 
09 


2 


781° 


17+0 
16-0 


14 

14 





7212 


03+0 
03-0 


04 
04 


^0 


Ds 





28lg 


02+0 
02-0 


02 
02 





681° 


01+0 
01-0 


04 
04 


1 


721° 


15+0 
14-0 


12 

12 





2712 


05+0 
05-0 


02 
02 


yB, 


D* 





231° 


02+0 
02-0 


03 
02 


1 


361° 


07+0 
07-0 


16 
14 


3 


631° 


23+0 
21-0 


21 

21 





9512 


04+0 
04-0 


06 
06 


At 


D' 





17lg 


01+0 
01-0 


01 
01 





94lg 


06+0 
05-0 


08 
08 


3 


5812 


23+0 
21-0 


22 
23 





8312 


06+0 
04-0 


06 
06 




D' 





14lg 


01+0 
01-0 


02 
01 





5llS 


01+0 
01-0 


04 
04 


2 


621° 


21+0 
19-0 


19 

19 





5312 


03+0 
03-0 


04 
04 


At 


D' 





121^ 


01+0 
01-0 


02 
02 





571° 


01+0 
02-0 


06 
06 


3 


1812 


19+0 
18-0 


18 
18 





6612 


03+0 
04-0 


04 
04 



pBcTIc 





61 


pBaVc 





61 


yBaJ/i' 





74 


ABaJ/i' 

^0 





53 


j^B^J/i, 





50 


j^B^J/i' 
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+0.03+0.01 
-0.04-0.01 
+0.03+0.01 
-0.04-0.01 
+0.01+0.03 
-0.01-0.03 
+0.01+0.02 
-0.01-0.02 
+0.01+0.02 
-0.02-0.02 
+0.02+0.02 
-0.03-0.02 



1.0912 
0.8612 
1.4512 



1.02 
0.80 



0.8112 



.00 + 0.05 
.02-0.05 
02+0.04 

03- 0.04 
03+0.09 

04- 0.08 
02+0.07 

02- 0.07 
00+0.05 
01-0.05 
02+0.05 

03- 0.04 



9912: 



1. 

1.18 



2.4612: 



2.39 
1.73 
2.22 



22+0.08 
20-0.08 
26+0.09 
24-0.09 
13+0.10 
13-0.10 
13+0.11 
13-0.11 
12+0.12 
12-0.12 
11+0.11 
10-0.11 



0.44 
0.17 
0.56 
0.50 
0.33 
0.45 



05+0.02 
05-0.02 
09+0.02 
09-0.02 
02+0.03 

.03-0.03 
02+0.02 
03-0.02 
01+0.02 
02-0.02 

.01+0.02 
01-0.02 



TABLE III: Results for the Be B, B* , Bs, B* form factors in the light front quark model. The uncertainties are from the Be 
decay constants and the decay constant of the final state mesons. 



F 



m 



rufit 



pB^B 

yB^B' 
aB^B' 
j\^BaB' 



n. (1.0+0.04+0.03 
'^•"■J-0.05-0.03 
p, f.q+0.04+0.03 
"•"■-'-0.05-0.03 
+0.17+0.32 



3.29 
0.47 
0.43 



n Qfi+O.05+0.08 
"■^"-0.07-0.07 
r, 01 +0.02+0.06 
"•oJ^-0.03-0.05 
+0.19+0.61 



1 1Q+0.09+0.01 
^•^="-0.09-0.01 
1 rr,+0.22+0.02 
^■0^-0.19-0.02 
+0.13+0.05 



+0.01+0.04 
-0.01-0.04 
+0.01+0.04 
-0.01-0.04 



0.68 
0.57 



+0.01+0.07 
-0.02-0.07 
+0.00+0.06 
■0.01-0.06 



2.65 
0.99 
1.16 



n. oq+0.04+0.01 
"■■^■J-0.04-0.01 

r, rr,+0.16+0.02 

"•■^^-0.10-0.02 
+0.27+0.10 



+0.04+0.04 
-0.04-0.04 
+0.07+0.03 
-0.07-0.03 



1.75 
0.31 
0.27 



+0.03+0.02 
-0.03-0.02 
+0.03+0.01 
■0.03-0.02 



pB,_B, 

yBcBi 

^0 

BaB' 



A, 



0.73 
0.73 
3.62 
0.56 
0.52 



+0.03+0.03 
0.04-0.03 
+0.03+0.03 
0.04-0.03 
+0.12+0.31 
0.15-0.29 
+0.00+0.04 
-0.01-0.04 
+0.00+0.04 



1.01 
0.87 
4.93 
0.75 
0.64 



+0.02+0.07 
0.04-0.06 
+0.00+0.05 
0.02-0.05 
+0.14+0.53 
0.19-0.47 
+0.00+0.07 
-0.01-0.07 
+0.00+0.06 



1.35 
1.77 
2.94 
1.13 
1.33 



+0.07+0.01 
0.08-0.01 
+0.24+0.04 
0.20-0.04 
+0.11+0.04 
-0.11-0.05 
+0.03+0.04 



+0.07+0.03 



0.35 
0.60 
1.78 
0.33 
0.28 



+0.04+0.00 
0.04-0.01 
+0.23+0.04 
0.14-0.04 
+0.25+0.07 
-0.21-0.08 
+0.03+0.02 



+0.03+0.01 



vanish, uncertainties of the continuum and the higher resonance interpolated by both of the axial- vector current and 
vector current are expected to be larger. In Ref. IJ], the authors also adopted the three-point QCD sum rules but 
different correlation functions are chosen. The form factors A2 " and A2 " " in Ref. have different signs with the 
other results. The large differences in different models can be used to distinguish them in the future. 

At the quark level, the P{V)lD decays are described as 6 ^ c{u)W~ c{u)W or c — > d{s)W^ d{s)l'^v. 

Integrating out the highly offshell intermediate degrees of freedom at tree level, the effective electroweak Hamiltonian 
for b ^ ulvi transition, as an example, is 

Heffih ^ ulvi) = ^K6S7^(1 - lz)bTr{l - l^W (34) 
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TABLE IV: Be — > D,D* and Be Ds,D* form factors at = evaluated in the literature. 







aB,d' 








DWf2]" 


0.154 


0.156 


0.145 


0.134 


0.224 


CNP[3] 


0.13 


0.05 


0.11 


0.17 


0.25 


NW[7] 


0.1446 


0.094 


0.100 


0.105 


0.175 


IKSfS] 


0.69 


0.47 


0.56 


0.64 


0.98 


Kiselevfg]' 


0.32[0.29] 


0.35(0.37] 0.43[0.43] 0.51(0.50] 1.66(1.74] 


EFGflO] 


0.14 


0.14 


0.17 


0.19 


0.18 


HZ[14] 


0.35 


0.05 


0.32 


0.57 


0.57 


DSVfie] 


0.075 


0.081 


0.095 


0.11 


0.16 




pBaD^ _ pBaD, 
^1 — ^0 


^0 








Kiselev[9]'' 
DSV[16] 


0.45[0.43] 
0.15 


0.47(0.52] 0.56(0.56] 0.65(0.60] 2.02(2.27] 
0.16 0.18 0.20 0.29 



"We quote the results with lu = 0.6 GeV. 

' The results out (in) the brackets are evaluated in sum rules (potential model). 



Since radiative corrections due to strong interactions only happen between the b quark and the u quark, they charac- 
terize the interactions at the low energy and the Wilson coefficient which contains the physics above the rrib scale is 
not altered. With the masses of leptons taken into account, the differential decay widths of Be Plv and Be VlD 
{I — e, /i, t) are given by: 



dT{Be PW) 



dTLiBc^VW) 



= ( 



VW) 



dq^ 



q 



1 



384m|;^7r3 



(35) 



= ( 



9 9 

1 - mf 2 



V X - {3mfX{ml^,ml.,q^)Al{q^ + 



384to| 7r3 



1 



2mi 



{m%^ - my - q^){mB, + my)Ai(q^) 



A2{q' 



'i \2 



X{ml^,m'y,q^)Gj,\Vc 



KM I 



384m|^7r3 



X < 



{mf + 2q^)\{ml^,m^y,q^) 



V{q') 



ms^ + my 



(m_B^ + my)Ai{q^ 



^X{m%^,m\,,q^) 



2\ 



(37) 



where the subscript +(— ) denotes the right-handed (left-handed) states of vector mesons. X{m^ ,m^,q'^) 



■>2N2 



Amg m^ with i = P,V. The combined transverse and total differential decay widths are given by: 



dq'^ dq^ 



dv- 



dV 
d^ 



dVL 
dq^ 



dVj 



(38) 



As we have mentioned in the above, the form factors A^"^ and ^4^""° only give small contributions to semileptonic 
Be decays. In these two channels, the two small variables m^ — my and q^ satisfy the inequality: q^ < g^^x = 
(ms^ — my)"^ <^ (™-Bc ~ fnv){mBe + my) = — m^. One can expand the decay width in terms of small variables. 
The variable A (m^^, my, g2) can be expanded as: X{m'^^,my, q"^) = {m\ — myY — 4,{m\+my)q^ + q'^ ~ (m'^^^ — my)^. 
From Eq. (|36p. we can see that the contribution from A2 to the longitudinal differential decay width contains a 



BaB' 
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TABLE V; Be — > r^c, J/tp, B, B* , Bs, B* form factors at = evaluated in the literature. 





1 


AB^J/ii 



AB^J/i, 
1 


.B,J/i, 
2 




DWf2]'" 


0.420 


0.408 


0.416 


0.431 


0.591 


CNPfS] 


0.20 


0.26 


0.27 


0.28 


0.38 


KTf4] 


0.23 


0.21 


0.21 


0.23 


0.33 


KLOfe]' 


0.66 


0.60 


0.63 


0.69 


1.03 


NWf7] 


0.5359 


0.532 


0.524 


0.509 


0.736 


IKS El 


0.76 


0.69 


0.68 


0.66 


0.96 


Kiselevfg]" 


0.66[0.7] 


0.60[0.66] 0.63(0.66] 0.69[0.66] 1.03(0.94] 


EFGflO] 


0.47 


0.40 


0.50 


0.73 


0.49 


IKS2^ 


0.61 


0.57 


0.56 


0.54 


0.83 


HNV[13l 


0.49 


0.45 


0.49 


0.56 


0.61 


HZ[14] 


0.87 


0.27 


0.75 


1.69 


1.69 


SDYflS] 


0.87 


0.27 


0.75 


1.69 


1.69 


DSV[16] 


0.58 


0.58 


0.63 


0.74 


0.91 




pB^B ptB^B 

1 


aB^b' 

^0 


^B^B- 




yB^B' 


DWf2]" 


0.662 


0.682 


0.729 


1.240 


5.690 


CNPfS] 


0.3 


0.35 


0.34 


0.23 


1.97 


NWf7] 


0.4504 


0.269 


0.291 


0.538 


1.94 


IKS [81'' 


0.58 


0.35 


0.27 


-0.60 


3.27 


Kiselev[9]'= 


1.27[1.38] 


0.55(0.51] 0.84(0.81] 4.06(4.18] 15.7(15.9] 


EFGflO] 


0.39 


0.20 


0.42 


2.89 


3.94 


AS[12] 




0.28 


0.17 


-1.10 


0.09 


HNVflS]" 


0.39 


0.34 


0.38 


0.80 


1.69 


HZ[14] 


0.90 


0.27 


0.90 


7.9 


7.9 


DSVfie] 


0.41 


0.42 


0.63 


2.74 


4.77 




pB,B, ^ pB^B, 
1 u 










DWf2]" 


0.715 


0.734 


0.821 


1.909 


5.657 


CNPfS] 


0.30 


0.39 


0.38 


0.35 


2.11 


CKMf5]^' 


0.403 


0.433 


0.487 


1.155 


3.367 


NWf7] 


0.5917 


0.445 


0.471 


0.787 


2.81 


IKSf8l'' 


0.61 


0.39 


0.33 


-0.40 


3.25 


Kiselevfg]'^ 


1.3[1.1] 


0.56(0.47] 


.69(0.70] 2.34(3.51] 12.9(12.9] 


EFGflO] 


0.50 


0.35 


0.49 


2.19 


3.44 


HNVfl3]'= 


0.58 


0.52 


0.55 


0.98 


2.29 


HZfl4] 


1.02 


0.36 


1.01 


9.04 


9.04 


DSVfie] 


0.55 


0.57 


0.79 


3.24 


5.19 



"We quote the results with u) = 0.6 GeV. 

' We quote the values where the Coulomb corrections are taken into account. 

■^The results out (in) the brackets are evaluated in sum rules (potential model). 

■^We add a minus sign to the form factors _Fi, Aq, Ai, A2 

'We add a minus sign for their predictions on the form factors. 

■^We quote the results which correspond to mj, = 4.9 GeV and uj = 0.4 GeV. 



12 



TABLE VI: Branching ratios (in units of %) and polarizations ^ of Be ^ Mlu decays. The first kind of uncertainties is from 
the Be decay constants and the decay constant of the final state mesons, while the second one is from the quark masses. The 
last kind of uncertainties is from the decay width of Be and the CKM matrix element Vub- The mass difference between an 
electron and a muon does not provide sizable effects in Be D^'hu and Be — > ric{J/ij)lv decays, but it does in Be B^*Hv 
and Be bI*Hv decays 





Be DeVe 


Be Dfll/fj 


Be DrVr 


Be rjeeUe 


Be ~* riefJ.U^ 


O /v 


n nn'?n+f ■0005+0.0005+0.0007 

'-'•""'^^-0.0004-0.0004-0.0007 


n nn'?n+0'Ooo5+o. 0005+0. 0007 

U.UUOU_oQoo4_o. 0004- 0.0007 


n nn9i +0.0003+0.0003+0.0005 

U.UUZl_o 0003-0.0003-0.0005 


n «7+0. 04+0. 04+0. 10 
'-'•" ' -0.07-0.04-0.10 


n «7+0. 04+0. 04+0. 10 
' -0.07-0.04-0.10 




Be rjcTVr 


Be — > BeUe 


Be BjJ.Pf_, 


Be BseUe 


Be Bs^U^ 


BTZ 


n 1 nn+0. 005+0. 014+0. 029 
U.i»U_Q 012-0. 013-0. 029 


n inQ+0-014+0. 013+0. 017 
U.1U9_Q 016-0.012-0.017 


n -1 n/i+0-013+0. 013+0. 016 

U.lU1_ooi5_o. 012-0. 016 


1 /lQ+0. 10+0. 15+0. 23 
^^•^^-0. 13-0. 14-0. 23 


-1 A -| +0.09+0.14+0.21 
-'-•^-'--0. 12-0. 14-0. 21 




Be D*eUe 


Be D*jj.i>^ 


Be D*TVr 


Be J/lpeUe 


Be J/ ^PflV^ 


BTZ 


n nn/ic;+0'0005+0. 0010+0. 0011 

U.UU'iO_o Q0Q4_Q 0008-0.0010 


n nn/it;+0'Ooo5+o. 0010+0. 0011 

U.UU4:0_o. 0004-0. 0008-0. 0010 


n nn97+o-ooo3+o. 0006+0. 0007 

U.UUZ ( -0.0002-0.0005-0. 0006 


■1 ^Q+0. 01+0. 15+0. 23 
^^•^^-0. 03-0. 14-0. 23 


■1 ^Q+0. 01+0. 15+0. 23 
1 •'±='-0.03-0.14-0. 23 




n co+0. 02+0. 02+0. 00 
U.DO_o. 02-0. 02-0. 00 


n eo+0. 02+0. 02+0. 00 
'^•'-'°-0. 02-0. 02-0. 00 


PI 7^1+0.01+0.02+0.00 
'J- "-'-0.01-0.02-0.00 


1 n/1+0-00+0.02+0.00 

^^•"^-0. 00-0. 02-0. 00 


1 n/1+0-00+0.02+0.00 
-■^•"^-0. 00-0. 02-0. 00 




Be —> J/lpTVr 


Be — » B'ei/e 


Be B*flU^ 


Be B*ei'e 


Be Bl^Ufj, 


BTZ 


n 070+0.002+0.042+0.056 
"••J "-'-0.005-0.038-0.056 


n -1^1+0.002+0.029+0.021 
U.iy:i_o. 004-0. 026-0. 021 


n -10/1+0.002+0.028+0.020 
U.10y:_Q 004-0.025-0.020 


1 Q«+0. 00+0. 34+0. 30 
l.aO_o. 03-0. 32-0. 30 


1 OQ+0. 00+0. 32+0. 28 
1.00_o. 03-0. 30-0. 28 




r, gi +0.01+0.01+0.00 
'-'•"-'--0.01-0.01-0.00 


1 r,7+0. 01+0. 02+0. 00 
^^•'J' -0.01-0.03-0.00 


1 nR+0-01+0.02+0.00 
-'-•'-"-'-0.01-0.02-0.00 


■1 -1/1+0.01+0.02+0.00 
-■-•-■-^-0. 01-0. 02-0. 00 


1 -1 -1 +0.01+0.01+0.00 
±•±±-0.01-0.02-0.00 



factor of \{m?g ,my,q^) while the Ai term is of the order yA(m^ ,my,q'^). Numerical results show that the ratio 

■7 ; fi-J — §^ — 2 TV and -, : — 3^ 5 jr is Smaller than 0.083 and 0.075 in the full region for q , 

respectively. It implies that the form factors ^^''^ and can be safely neglected in the decay width ^. 

Integrating the differential decay widths over the variable q^, one obtains partial decay widths and polarization 
fractions. The lifetime of the Br meson and the relevant CKM matrix elements are used as 



2M- 



TB, - (0.46 ± 0.07)ps, \Vcb\ = 41.2 x 10-^ = (3.93 ± 0.36) x 10"^, = 0.230, \Vcs\ = 0.973, (39) 

where the small uncertainties in the other CKM matrix elements are neglected. Our predictions of branching ratios 
and polarization quantities in semileptonic decays are given in table IVll The three kinds of uncertainties are 
from: the decay constants of the B^ meson and the meson in the final state; the constituent quark masses; the lifetime 
of Be together with the CKM matrix elements. The first kind of uncertainties in the Be {D, Dg, J /'ip)li> decays is 
very small, as the uncertainties in decay constants of D and J /'ip are small. The different mass between the electron 
and muon does not have sizable effects on — )■ u, c semileptonic decays, but the branching ratios of c — > u, s transitions 
are altered by roughly 5%. Branching ratios of Be Plv decays are smaller than the corresponding Be Vlv ones, 
partly because there are three kinds of polarizations for vector mesons. Among the four kinds of transitions at the 
quark level, there is an inequation in chain: 

BTliBe D*lv) < BTZ{Be B*lv) < BTZ{Be J/^lv) < BTZ{Be B*Jiy), (40) 

where we have taken decays involving a vector meson as an example. To understand this inequation, three points are 
essential. The CKM matrix elements for these four kinds of decays are given as: 

\Vub\<t:Veb<t:\Ved\<t:Ves. (41) 

The form factors at zero- recoiling point roughly respect: 

F{Be ^ D*) < F{Be ^ J/ij) ^ F{Be ^ B*) ~ F{Be ^ Bl). (42) 



^ In Be — > {B* ,B*)V decays, the analysis is similar: is replaced by the mass square of the vector meson rUy 
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The phase spaces in Be D* and Be J/tp transitions are much larger than those in Be B*,B* transitions, 
which can compensate for the small CKM matrix element in Be J/tplD decay. These predictions will be tested at 
the ongoing and forthcoming hadron colliders. 



IV. CONCLUSION 



Due to the rich data, measurements on the CKM matrix elements are becoming more and more accurate. Be meson 
decays provide another promising place to continue the errand in B meson decays. They also offer a new window 
to explore the structure of weak interactions. Although the Be meson can not be produced on the two B factories, 
it has a promising prospect on the ongoing and forthcoming hadron colliders. Because of these interesting features, 
we have studied the Be transition form factors in the covariant light-front quark model, which are relevant for the 
semileptonic Be decays. 

Comparing our predictions with results for the form factors in the literature, we find large discrepancies which 
may be useful to distinguish various theoretical methods. Our results for the form factors A2 in Be — » B* and 
Be B* transitions strongly depend on the decay constants of the B* and B* mesons, which gives large theoretical 
uncertainties to the form factors. For Be BP decays, the relevant form factor Ag is almost independent of A2: 
Aq ~ Ai. For semileptonic Be decays (also Be B*V decays), contributions from A2 are at least suppressed by 
a factor of 0.08 compared with those from Ai. Thus the large uncertainties from A2 will not affect the physical 
observables. 

Be D, D* , Ds, D* form factors at maximally recoiling point are smaller than Be ije, J/tp, B, B*, Bs, B*, while 
the Be D, Ds, rje form factors at zero recoiling point are close to each other. The SU(3) symmetry breaking effects 
in Be D, Ds and Be D* , D* are quite large; but in Be — > B, Bg and Be ^ B* , B* transitions, the SU(3) breaking 
effects are not large. Semileptonic Be — > {rje, J /ip)!!^ and Be {Bg, Bl)lv decays have much larger branching fractions 
than the other two kinds of semileptonic Be decays. In the three kinds of Be — > Vlv decays, contributions from the 
longitudinal polarized vector is comparable with those from the transversely polarized vector. These predictions will 
be tested at the ongoing and forthcoming hadron colliders. 
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APPENDIX A: RELATIONS OF DIFFERENT DEFINITIONS OF FORM FACTORS 

In the literature, various conventions for the Be ^ V form factors have been adopted. In this appendix, we will 
collect their conventions and compare them with the BSW form factors. In Refs. 0, Q, 0, 0] , the authors defined the 
Be^V form factors as: 

{V{P",e")\V^\Be{P')) = -e^.o.pe"*''P^q''Fv{q% (AI) 
{V{P\e")\A^\Be{P')) = iF^iq^y;: + iF+{q^){e"* ■ P)P^+ iF^{q^){e"* ■ P)q^, (A2) 

These form factors are related to the BSW form factors by: 

= {mB^+mv)Fv, = ^ , A^"" = -{ruB^ + mv)F+, (A3) 
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^ m^^ + my ^py ms. - my ^^^^ 



The definition of form factors g, f,a-^-,a^ in Ref. (7| is similar with ours in Eqs.([T][3|) except for a phase i. In 
Ref. Qlli, the following definition for the form factors is adopted: 

{V{P",e")\V^-A^\B,{P')) = eT{'9'"'P-qA„ + P^'P''A++q''P''A-+te^'''P''Ppq^V), (A5) 

niB^ + my 

where A+ corresponds to the BSW form factor and their form factor A^^^^ is related to the BSW form factor 
Af^: 

Ar = ^""''^'""--"^^^ (A6) 
niB, + my 

In Ref. [3], the Be V form factors are defined as: 



{V{P'\ e")\V^ - A^\B,{P')) - -te'l*{mB^ + my)A, + iP^{e"* ■ q) f + tq^{e"* ■ q) 



mB^+my ms^+ray 
niB, + my 

The form factors Af^ and V^^ are the same with the relevant BSW form factors; their form factor corresponds 
to the BSW form factor A^^ . 



APPENDIX B: SOME SPECIFIC RULES UNDER THE p" INTEGRATION 



When performing the p~ integration, one needs to include the zero-mode contribution. This amounts to performing 
the integration in a proper way in this approach. To be more specific, for p[ under integration we use the following 
rules 



22, 



23| 



= P^A['^ +q^A^^'>,N2~^ Z^, 



(1) 



l(2) 



(2) 



l(2) 



(Bl) 



where the symbol = reminds us that the above equations are true only after i nteg ration. A^'-* are functions of xi,2, 
p']^, p'^ ■ q± and q^, and their explicit expressions have been studied in Ref. 22l |23||: 



N[ + -ml + {l^ 2xi)M'^ + {q^ + q ■ P) 



,(1) _ xi (1) _ (1) p1 ■ q± (2) _ ,2 {p'±-q±? 

^1 - V 2 -^1 "2 '"^1 ^ P-i- 72 ■ 

z y q 



aC^) _ A^'-> A'^'> /l^^-* - (A 



(1)^(1) /l(2) 



q 



(B2) 



We do not show the spurious contributions in Eq. (|B2p since they are numerically vanishing 
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